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ABSTRACT 



We present a semi-analytical model in which galaxy collisions and strong 
tidal interactions, both in the field and during the collapse phase of groups and 
clusters help determine galaxy morphology. Prom a semi-analytical analysis 
based on simulation results of tidal coUisions (Aguilar & White 1985), we 
propose simple rules for energy exchanges during collisions that allow to 
discriminate between different Hubble types: efficient collisions result in the 
disruption of disks and substantial star formation, leading to the formation 
of elliptical galaxies; inefficient collisions allow a large gas reservoir to survive 
and form disks. Assuming that galaxy formation proceeds in a = 1 Cold 
Dark Matter universe, the model both reproduces a number of observations and 
makes predictions, among which are the redshifts of formation of the different 
Hubble types in the field. When the model is normalized to the present day 
abundance of X-ray clusters, the amount of energy exchange needed to produce 
elliptical galaxies in the field implies that they formed at z ^ 2.5 while spiral 
galaxies formed sX z < 1.5. The model also offers a natural explanation for 
biasing between the various morphological types. We find that the present 
day morphology-density relation in the field is well reproduced under the 
collision hypothesis. Finally, predictions of the evolution of the various galaxy 
populations with redshift are made, in the field as well as in clusters. 



Subject headings: cosmology: theory - galaxies: formation - galaxies: 
morphologies - galaxies: kinematics and dynamics 
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1. Introduction 

Gravitational interactions between galaxies are believed to play a major role in 
determining galaxy physical properties (Spitzer & Baade 1951; Toomre & Toomre 1972; 
Toomre 1974; Schweizer 1996; Schweizer & Seitzer 1992; Cole et al. 1994). For example, 
galactic interactions are likely to enhance star formation rates and could explain the intense 
star formation seen in IRAS galaxies. Indeed, these highly luminous infrared galaxies often 
show evidence of undergoing tidal interactions or mergers (Clements et al. 1996; Leech et 
al. 1994). 

On theoretical grounds, models have shown that close encounters of galaxies stimulate 
cloud growth and star formation. Stellar and mass distributions in galaxies are likely to be 
altered by tidal collisions, and, as a result, the morphology of a galaxy may evolve. For 
example, the disks of spiral galaxies might be puffed up by such tidal encounters, and it has 
been suggested (Richstone 1976) that spiral galaxies might be converted into SO or elliptical 
galaxies by such a mechanism. Following this idea, galaxy morphologies might be largely 
determined by the local properties of the environment in which they form and evolve, since 
tidal collisions between galaxies occur more frequently in denser regions of the universe. 

Alternatively, galaxy morphologies may be related to intrinsic properties of the 
primordial fluctuation field. In the Cold Dark Matter (CDM) theory in which the density 
parameter Qq — 1, a bias must be invoked whereby the overdense regions contain galaxies, 
and the underdense regions (voids) are deficient in luminous galaxies, so that the luminous 
matter has stronger density correlations than the underlying dark matter. Early discussions 
of CDM used Gaussian fluctuations as described by the peaks formalism (Bardeen et al. 
1986) and introduced a bias of the galaxy types with respect to mass by identifying 3a 
primordial fluctuations with protoelliptical galaxies and 2a fluctuations with protospirals 
in order to match the observed frequency and clustering of luminous galaxies (Blumenthal 
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et al. 1984; Evrard 1989; Evrard, Silk & Szalay 1990). Though this approach offered the 
advantage of reproducing several observed correlations, no physical mechanism was proposed 
to explain this intrinsic biasing until Dekel & Silk (1986) argued that la fluctuations that 
satisfy the cooling criterion for galaxy formation (Rees & Ostriker 1977; Silk 1977) would 
be systematically of low-mass and shallow gravitational potential wells so that the rms 
density fluctuations should be especially vulnerable to disruption by supernovae explosions. 
A prediction of this model is that an extensive distribution of 'failed galaxies', identified as 
dwarf elliptical galaxies, would populate the low-density regions of the universe. However, 
it has been argued that observations strongly constrain the hypothesis that dwarf ellipticals 
can account for biasing (Bingelli 1989), the more luminous, observed dwarf ellipticals 
clustering together with the bright galaxies. 

An alternative physical bias mechanism may naturally involve galaxy interactions as 
far as these interactions are assumed to trigger star formation. According to this approach, 
low-density regions of the universe should contain nascent or unborn galaxies, essentially 
gas clouds, since collisions are much less frequent in underdense regions. Conversely, in 
denser regions, substantial interactions at the epoch of formation would induce early star 
formation and subsequent collisions may redistribute the galactic stellar content and alter 
the galaxy morphology. 

In this paper, we build a simple semi-analytical model of galaxy interactions to 
understand galaxy morphologies and how they relate to the properties of the environment 
and the fluctuation field at the epoch of formation. We shall show that simple collision 
rules provide an explanation of the morphological distribution of galaxies with respect to 
environment, as well as of fundamental correlations in their structural properties. Moreover, 
biasing of galaxies of different morphological types with respect to mass originates naturally 
in our scenario. The model is based on the results of N-body simulations of tidal collisions 
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(Aguilar & White 1985; AG85 hereafter). The paper is organized as follows: we first 
describe our model. Galaxy collisions are characterized in terms of the rate of change of 
binding energy induced in a given galaxy, at a given epoch, by tidal encounters with a set 
of background galaxies (§2). The cases for field and cluster galaxies are treated separately 
as the physics of collisions is quite different in these two environments. In §3, the total 
change of binding energy in a galaxy that went through a series of tidal collisions since 
its formation (we shall refer to this quantity as the 'collision factor') is computed and 
characterizes the coUision history of the galaxy. Scalings of the coUision factor with local 
galaxy density are obtained. In §4, wc propose a phenomenological definition of galaxy 
morphological types that can be easily expressed as conditions on the collision factor (§4.2). 
We show (§4.3) that these conditions are in turn conditions on the formation redshift of 
galaxies, and we define redshift cuts, in the space of formation redshifts, that delineate 
spiral, SO and elliptical galaxies. Sections 5 and 6 are devoted to normalizing our model to 
various sets of observations, assuming CDM initial conditions. Predictions of the model are 
given in the remaining sections. The model predicts the redshift of formation (§7) as well 
as the relative bias between galaxies of various morphological types (§8) in the field. The 
morphology- density relation is well reproduced by the model (§9). The predictive power 
of the model is illustrated in §10, where the evolution of the morphological populations 
with redshift is predicted, under CDM initial conditions, in the field as well as in denser 
environments. Throughout the paper, the case for an Einstein-de Sitter [Qq = 1) with a 
Hubble constant Hq — 50 km/s/Mpc is assumed. 

2. Characterization of galcixy collisions 

We characterize galaxy collisions in terms of the rate of change of binding energy E of 
a galaxy colliding with an other. This can be expressed in a general form for a test galaxy 



- 6 - 



(of mass M and radius R) interacting with a set of background galaxies (of mass Mp and 
radius Rp) with relative velocity v as (Richstone 1975; AG85): 

K = npvR^fE{Mp/M,V/v) (1) 

where A stands for d\nE/dt. Hp is the number density of background galaxies, and V the 
internal velocity dispersion of the test galaxy. The dimensionless function Je is plotted in 
figure 8 of AG85 for equal masses galaxies. It is easy to see from this figure that scales 
roughly as the square of the ratio V/v so that equation (|^) simply reads as: 

K = fR^{^)\np (2) 

where / is the slope of the function fE and is equal to / ~ 30 (AG85; figure 8). For unequal 
galaxy masses, two modifications of @) may be anticipated: 1) R should be replaced by 
i?>, where i?> = max(i?, and 2) powers of the mass ratio M/Mp may appear. For a 
test galaxy (_R, M) perturbed by a set of background galaxies (-Rp, Mp) with distribution 
dup = f]{Mp, z)dMp, the perturbation 6V on the stellar velocity field can be estimated from 
the impulse approximation (Spitzer 1958) and straight line trajectories (provided v ^ V) 
to be: 

SV GMp R 

TT T7 — 2 ^-^f 

V vV p"^ 

where p is the impact parameter which for the dominant contribution to the rate of energy 
exchange can be taken to approximately equal i?>. Using the virial theorem for the test 
galaxy, equation (^ becomes: 

6V V R^ Mp 

T ^ ~W,~M ^ ^ 

For equal mass galaxies, we have from (^ that 5V/V ~ V^/f so that A scales as {5V/VY - 

equation (0). We assume this still holds for unequal masses with 5V/V given by (^. The 

rate for energy exchange now reads as: 
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where the integral is performed over the masses of background galaxies. As equation (j^) 
relies on the impulse approximation and straight line trajectories, it is valid only for v ^ V. 
This latter condition is fulfilled in galaxy clusters, for example, as the relative velocity 
between galaxies is expected to be of the order of the cluster velocity dispersion a ~ 1000 
km/s, whereas the internal velocity dispersion of galaxies is of the order of ~ 200 — 300 
km/s. In the field, the relative velocity of colliding galaxies may be of the same order as 
their internal velocity dispersion and the relative trajectories can not be approximated 
by straight lines. The focusing due to mutual attraction can be easily incorporated into 
equation (^). If p and v are the initial impact parameter and relative velocity, energy 
conservation during the collision implies that, in the frame of the reduced particle of mass 
/i = MpM/{Mp + M): 

1 2 1 2 GMMp 

-fiv = -fiv^^i (6) 

^ ^ J^col 

where Vcoi is the relative velocity at closest approach Rcoi, and G is the gravitational 
constant. Conservation of momentum per unit mass L implies L = pv = RcoiVcoi from which 
we get: 

, _ „2 M , 2G(M + M,) ^ „2 . , 2G(M>) 

Using the virial theorem for the largest of the two colliding galaxies (of mass M> and 
internal velocity dispersion V^>), we have: 



P ~ ^col 



V 



Taking R^oi ~ -R>, gravitational focusing roughly amounts to multiplying Ri in equation 



by the factor 

i + ( 

For a similar reason, the relative velocity entering equation (^) is increased by the square 
root of (^. However, the energy taken out of the relative motion to be injected into the 
internal stellar motions lowers the relative motion in a similar way. If f ~ V", the two effects 
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nearly cancel each other. We do not attempt to model these fairly complex processes, but 
simply use the rate modified by the factor (P). The scaling with mass of the latter is 
needed in order to take into account the strong attraction produced by large galaxies, and 
will be important for the scaling of the rate of change of energy with mass and redshift. 

We now integrate (|^) over the masses of the background galaxies. We use i?> = R, the 
radius of the test galaxy, so that the test galaxy is larger, on average, than the background 
galaxies. We define / Mprj{Mp)dMp = ep, where p is the mean density in the universe and e 
is the mass fraction in galaxies, so that we get: 

and 

where M in equation ([10|) and (|ll]) is defined by: 

J Mj,ri{Mp)dMp- ^ ' 

M is dominated by the larger masses due to the Mp factors: hence, we expect M to be close 
to the mass of the larger galaxies. 

To summarize this section, we have obtained scaling relations for the rate of change of 
binding energy as a function of the various parameters characterizing the colliding galaxies. 
Table 1 gives a list of the names and meaning of the parameters used so far. Table 2 
summarizes the results we shall build our model on in the next sections, and the various 
assumptions used. 
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3. Collision history of galaxies 

The rate of change of binding energy for a galaxy interacting with surrounding galaxies, 
between redshift z and 2; + is 

• dt 

dlnE = A—dz, (13) 
dz 

where A has been given by (|T0|) and (jl^) in the previous section. Integrating (|1^) over time 
from the redshift of formation Zni of the galaxy up to redshift z yields the total increase of 
binding energy due to a series of collisions experienced by the galaxy during its lifetime. 
This integrated quantity characterizes what we shall call the 'collision history' of a galaxy. 
To evaluate this quantity, we need to use the relevant scaling with redshift of the various 
quantities entering equations ([T0|), (|ll|) and (p!2D , which have been listed in Table 1. These 



scalings will be different in cluster and field environments. We shall therefore distinguish 
from now on the cases for cluster and field galaxies. We remind the reader that all the 
results given in this paper are for an Einstein-de Sitter universe, that is Qq = 1, and for 
h = 0.5, where h is the Hubble constant in units of 100 km/s/Mpc. 



3.1. Redshift evolution of M 

By construction (equation (0)), M depends on the mass distribution ri{M) of the 
background galaxies. rj{M) evolves with redshift and so does the average mass M in our 
model. We evaluate the scaling of M with redshift in this section. 



3.1.1. Press & Schechter mass function 

With the assumption that the initial fiuctuations are Gaussian distributed, the mass 
function of dark halos can be obtained from the condition that a fiuctuation is non-linear 
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at a given mass scale, but not at an immediately larger scale (Press & Schechter 1974; 
Schaeffer & Silk 1985). The number density of non-linear condensations of mass between 
M and M + dM at redshift z is, in an Einstein-de Sitter universe: 

where po is the present day average density of the universe and 5c is the linearly extrapolated 
threshold on the density contrast required for structure formation. We adopt the canonical 
value of the spherical model (Lemaitre 1933; Peebles 1980; Gunn & Gott 1972), i.e. 
5c ~ 1.68, and the average density of an object collapsing at redshift z is 178po(l + zY . The 
relative mass fluctuation 5M/M in a volume that contains a mass M in the linear stage 
enters equation (0) through its variance 

^WH(f)?. (15) 

which is known as soon as the primordial fluctuation spectrum is specified. In the following, 
we shall use the Cold Dark Matter (CDM) spectrum that can be parameterized (e.g., 
Narayan & White 1988) as a function of the comoving scale i?, corresponding to the mass 
scale M, as 

(Jcdm{M) = 16.3(1 - 0.3909i?°-^ + 0.4815i?°-^)"^7^ (16) 

where the bias parameter b has been introduced and is specified by the amplitude of 
underlying matter fluctuations at 8h~^ Mpc, erg- The parameterization ([T6|) is given for 
h = 0.5, a value that we shall adopt throughout the paper. 

At a given redshift z, we shall now require that gas cooling occurs after virialization to 
allow for gas fragmentation and star formation (Silk 1977, Rees & Ostriker 1977). This is 
important for evaluating the number density of background galaxies, as only a subset of 
all mass condensations counted for by the Press & Schechter prescription ( |T^ ) are actual 
galaxies. 
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3.1.2. The cooling constraint 

The condition for cooling to be effective is a condition on the relative importance of 
the cooling time-scale and the dynamical time-scale. Various processes contribute to gas 
cooling. Assuming that during the collapse, the gas is shock heated and reaches virial 
temperature before settling into a disk, we are led to consider only those cooling processes 
that are efficient at temperatures of the order of the virial temperature. For simplicity, we 
assume that line cooling is the dominant mechanism, and shall adopt the following cooling 
function for gas temperatures of the order ~ 10^ — 10^ K (e.g., Sutherland & Dopita 1993): 

A(r) ^ 2.5 X lO^^^T^^/^n^ ergcm'^s"^ (17) 

where n is the particle density of the gas within a galaxy and T the virial temperature. 
Zero metallicity has been assumed. The typical cooling time-scale is then 

where k the Boltzmann's constant. If the gas makes up a constant fraction Tb of the total 
mass M of the galaxy, and is uniformly distributed within the virial radius R, then the gas 
density n is 

3 J^bM _ _ ,/ M N/ R \-^rJ^E 



n = ^ = 1.6 X lO^^f— )( — ) ( — ) cm"^ 19 

47r/imp R^ V1012Mq/ UOOkpc/ VQ.l/ ^ ^ 



where is the proton mass and /i the mean molecular weight, that for an hydrogen- helium 
plasma with primordial abundances is /i ~ 0.6. The temperature of the gas is obtained 
from the viral equation: 

kT ^ limp— ^ ^^'^p~^ (20) 



where V is the 3-D galaxy velocity dispersion. In (|20|) , we have assumed that the galaxy is 
spherical and the gas is homogeneously distributed within R. Plugging (|l^) and (|3) into 
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we obtain the following cooling time-scale: 

tcooz - 1.7 X 10 (-) (t^^t^) (t^) yr (21) 



The dynamical time-scale can be estimated from the time taken for the galaxy to collapse 
after turn around: 



/ . r- .„q/ M \-l/2/ R x3/2 

tdvn^-\ — —^l.bxlO^i—- ) ( ) yr (22) 

where Rta ~ 2i? is the turn around radius, and G the gravitational constant. For efficient 
star formation, we require (|2l|) to be smaller than (^2]) implying: 



M < M, ^ 9 X 10^^ Mq (23) 

where is a critical mass whose value is fixed by the physics of cooling. 

For the purpose of evaluating the number density of background galaxies, we shall use the 
constraint ( p3|) that gives an upper limit of the mass of galaxies. In figure 1, we plot the 
average mass M, computed from ([l^), as a function of redshift (solid line) when the cooling 



constraint is applied. Integrals entering ( [I^ ) have been computed using a minimum mass 
Minf = 10^° Mq. For comparison, M obtained from (|I4|) with no cooling constraint is also 
plotted (dashed line). When the cooling constraint is used, the average mass M exhibits 
only a weak dependence on redshift over the range z ^ 3 — 4, much weaker than when no 
cooling is used. Moreover, the average mass is found to be of the order of M*. We shall 
adopt this value hereafter and take M constant with redshift for z < 3 — 4. 



3.2. Field galaxies 



The average density of the universe in equation (|TT]) scales with redshift as 

p = Po{l + zf 



(24) 
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where Pq is the average density of the universe at the present epoch. We assume that the 
average density in galaxies reflects the density of the universe at the epoch of formation 
Zni within a universal factor (that depends on the properties of the collapse). We then can 
relate the radius R of the test galaxy to its mass M as 

For simplicity, we have expressed the relevant quantities for the test galaxy in terms of the 
same quantities for the typical M^, galaxy introduced in the previous section. A typical M^, 
galaxy has radius i?*, mass M^,, internal velocity dispersion V^, and form at redshift z^. 
Using equation (^) and the virial theorem, we obtain the following scaling of the internal 
velocity dispersion of the test galaxy V with redshift: 



We finally need to know how the relative velocity v evolves with redshift. We assume that 
the mean relative velocity of colliding galaxies can be inferred from the peculiar motions of 
galaxies in the linear regime at redshift z: 

= = voil + z)-^/^ (27) 

where x is a comoving coordinate, a is the scale factor of the universe, and fo the present-day 
galaxy peculiar velocity in the field. This assumption is valid at high z in the linear regime. 
An accurate scaling at low redshift should take into account pairwise velocity on small scales 
in the non-linear regime; however, (p7D introduces only a weak dependence on redshift, and 



we shall adopt this scaling in the following. Plugging equations (p^) , (p5|) , (p6D and (|27|) into 
equation ([TT|) , we get 

A = efRfj^{^)\l + zf/^. (28) 

The redshift integration of equation (^) from the redshift of formation Zni of the test 
galaxy then yields the total increase of binding energy due to collisions, up to an epoch 
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characterized by redshift z. In the following, we shall call this quantity the 'collision factor' 
and denote it A. In an Einstein-de Sitter universe, redshift and time are related through: 



dt 



(29) 



where Hq is the present day value of the Hubble constant. Integration of (|28|) over redshift 
thus leads to: 



A^'^^'^iz) = I A^dz = AJ1 + z^if 



^dt 

ii dz 

where A* is the dimensionless constant 



1 



1 + -2 \3 



1 + Znl 



(30) 



A. = e(//3)i?^p) 



vj HoM,' 



(31) 



A numerical estimate of (|3T|) can be obtained by fixing the parameters describing a typical 
M^, galaxy today, and the relative velocity Vq. 



3.3. Cluster Galaxies 

For field galaxies, collisions are only important at early times when the density is large. 
Similarly, we expect little effect in the lower density regions of clusters. In the very high 
density regions of rich clusters, however, the energy exchange due to collisions is high and 
is given by equation ( p!0D with p replaced by the local density p and v by the local galaxy 
velocity dispersion of the cluster. Also, in the focusing factor (^, the term (Vy/vY) that 
was dominant in the field can now be neglected. Using equations (p5|) and (|2^) , equation 
(PUP thus becomes: 

Performing the same integration over redshift as in the previous subsection, now for cluster 
environment, yields the change of biding energy for a cluster galaxy from the epoch Zni up 
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to redshift z\ 



A 



with 



cluster 



-3/2 



(33) 



(34) 



where the average density of the universe, po, has been introduced for convenience. 



3.4. Scaling of A'^'"'^*'^'' with density 

In equation (P5|), both p and v can be related to the local galaxy number density n, 
provided a model is assumed for the density distribution of the cluster. Let us consider for 
instance the following cluster profile: 

p{R)=Pc{§y (35) 

where pc is the cluster density at a radius Rc from the center, and depends on the mass 
and redshift of formation of the cluster. By use of the virial theorem and (|35D, the velocity 
dispersion v can be expressed as a function of the cluster density p: 

- = -4-)H-)^ (36) 

where Vc is the cluster velocity dispersion at radius Rc- As expected for p = 2, the velocity 
dispersion of a given cluster is constant and independent of the density p. This velocity, 
however, is not the same for all clusters since it depends on p^. 

v = vJP^f' p = 2. (37) 

The collision factor (|33|) scales, in this simple 

^cluster ^ ^^-1/2 _ ^gg^ 
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We can relate the cluster density p to the local galaxy number density n by assuming that 
galaxies trace mass, that is n oc p. We shall however adopt a somewhat more realistic model 
for the cluster density profile, namely an isothermal (3 model (Cavaliere & Fusco-Femiano 
1976) in which the galaxy number density in the cluster is related to the gas density profile 
PgiR) by: 

Pg{R) oc n^(i?) (39) 

where /? is a fitting parameter to the projected X-ray surface brightness of the cluster 
(Jones & Forman 1984). Assuming that the gas makes up a fraction of the total mass that 
is roughly constant throughout the cluster, we have for the cluster density p: 

p{R) oc pg{R) oc n^{R) oc (^)~^^ R > Rc- (40) 



At large radius, p behaves with radius as equation (|35|) with p = 3/3. Expressed in terms of 



the galaxy number density, equation (pSD now scales, for collisions occurring at i? > R^ as: 

^cluster ^ ^/3^-/3/2 (4^) 

where Uc is the galaxy number density at radius Rc- We are led to two different cases: 1) 
collisions at a fixed density n within a given cluster, and 2) collisions at a fixed density 
averaged over a set of clusters. The scaling with density of the collision factor (^l]) will be 
quite different in these two cases. It is of interest of considering case 2) for the purpose of 
comparing the predictions of our model to observations, which we will present later in this 
paper. For example, (^) can be averaged, at a fixed galaxy number density n, over a set 
of clusters with varying 'central' density ric- Bahcall (1979) shows that the distribution of 
Abell clusters decreases as a function of luminosity L according to a Schechter law 

^(L) oc exp ( - L/La) (42) 

where La is the cluster optical luminosity within the Abell radius Ra- Assuming that the 
cluster optical luminosity is proportional to the galaxy number density, the central density 
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distribution of Abell clusters is 

oc n"^ exp ( - nj ria^ (43) 

with Ha oc La/R\- The colhsion factor (^Tj) averaged at fixed n over the distribution (|43|) 
for ric > "n, is then 

< A'^'"^*'^'' >~ < n;^/2 > (44) 
where <> denotes average over (^31): 



< n-/'^ >= n-^/'fin/n,) (45) 

where the function /(x) is defined by: 

_ r jrt'-'^/'exp(-i)A j 

' L /r*"exp(-()<i( i- * ' 

f{n/na) is dominated by the contribution of densities n na and behaves roughly as 
{n/na)~^^'^ ■ Thus, (|^) becomes 

< A^'""*^" >~ n^/l (47) 



In the following, we will use equation (^Tj) for the purpose of evaluating the morphology- 
density relation in our model. 



3.5. Comparison between the field and clusters 

Several differences in the collision factor in the two environments considered originate 
from the different scalings used in these two cases: 

• In the field, the collision factor depends essentially on the present day density of the 
universe, whereas in clusters, A is a function of the local density p. As a consequence, 
the energy exchange in clusters are generally much more important than in the field, 
and depend on the environment, since they increase with increasing density. 
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In the field, collisions at high redshift are more efficient (in terms of the fractional rate 
of change of binding energy per collision (pSD) than recent collisions. This is a direct 
consequence of the strong increase of the average density of the universe with redshift. 
This effect overseeds the fact that collisions for galaxies with small radii (that is, at a 
given mass, those formed at high Zni - equation (|25|) ), are less efficient than for large 
galaxies. Thus, at a given mass, ( pi]| ) implies that the higher the redshift of formation 
of a galaxy, the higher the collision factor. On the contrary, in a given cluster, the 
fractional rate of change of binding energy per collision (|32D is constant in time, as 
soon as the cluster has virialized. At a given mass, the higher the redshift of formation 



of the galaxy, the lower the galaxy radius (ESf), and the lower the collision factor (|3 



4. A phenomenological definition of Hubble types 

4.1. Angular momentum and morphological types 

A key phenomenon to understand when elaborating a consistent picture of galaxy 
formation is the angular momentum history of galaxies of various Hubble types. It is 
believed that the angular momentum of primeval galaxies originates from tidal torques 
induced by the presence of neighboring galaxies (Doroshkevich 1970; White 1984). At 
the early epochs of galaxy formation, protogalaxies indeed interact strongly with their 
surroundings. The angular momentum acquired at that time may have a definitive impact 
on galaxy morphologies. It would be however misleading to associate isolated structures, 
i.e., objects that are less subject to tidal torquing, with elliptical galaxies. Following 
this argument, one would then expect to find fewer low-spin galaxies in high density 
environments, a prediction in disagreement with the observation of an increasing fraction of 
early-types when one goes towards richer clusters (Dressier 1980). 
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It is unlikely that the large amount of spin needed to rotationally support fully grown 
disk galaxies may be explained solely by the primordial angular momentum (e.g., Barnes & 
Efstathiou 1987). Dissipative collapse of the radiatively cooled gas in a dark halo must be 
invoked. Assuming that the gas angular momentum is conserved during the collapse phase, 
one needs a collapse by a factor 10 to achieve rotational support (Fall & Efstathiou 1980). 
A large initial angular momentum is not a sufficient condition for the formation of disk 
galaxies, as various processes, such as tidal interactions occurring during the lifetime of a 
galaxy, may induce angular momentum transport from the galaxy center to the outer parts. 
Thus an initially rapidly rotating protogalaxy could end up as an elliptical galaxy. Recent 
simulations indeed find an extremely high efficiency of angular momentum transfer resulting 
in disks that are far too small. Even if the gas starts to sink towards the galactic center, 
the settlement into a disk structure requires a gentle infall that could be perturbed by the 
tidal field. Indeed, modeling of disk formation as well as chemical evolution models (Lacey 
& Fall 1983; Rocca-Volmerange & Schaeffer 1990) supports slow disk formation via gas 
infall from a pre-existing halo. Spirals are found to have a relatively constant rate of star 
formation over the past 10 Gyr, and infall provides a possible means of regulating the gas 
supply and maintaining the disk in a state of marginal gravitational instability (Sellwood 
& Carlberg 1984). However, as soon as a protogalactic cloud collapses and decouples from 
the universal expansion, collisions are expected to occur due to its relative velocity with 
respect to other clouds. Virialization should be effective within one or a few Hubble times 
after formation and collisions during this period are not expected to greatly modify the 
final structure. More recent collisions, however, would inhibit the gentle infall of the gas 
needed to allow for the formation of a disk. Clearly, the collision history over the whole life 
of a galaxy since its epoch of formation should be taken into account if one wants to predict 
its morphological type. 
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4.2. Definition of Hubble types 

The previous discussion leads us to define, in the framework of our model, the 
morphological types as follows: the condition for a spiral galaxy to form out of a cloud 
that first became non-linear at a redshift Zni is to experience few, if any, collisions between 
the epoch of its formation and the epoch under consideration. Strong collisions will, on 
the contrary, prevent the gas from settling into a disk and allow for tidal exchanges that 
average out angular momentum. We shall consequently assume that a sizeable number of 
collisions between Zni and the epoch characterized by z leads to the formation of elliptical 
galaxies. These rules for the formation of morphological types can be expressed in terms of 
conditions on the collision factor A, for field and/or cluster galaxies. At any redshift z, the 
collision factor for a spiral galaxy is assumed to satisfy the following condition: 

A(^)<AX (Spiral) (48) 

where A^^* is a redshift independent threshold on the collision factor and is a free parameter 
in this model. An elliptical galaxy will be such that: 

A(^) > Affl (Elliptical) (49) 

where A^" is another free parameter of our theory such that Affj: > A^^*. Finally, galaxies 
whose collision factors satisfy 

Ar<A(^)<A,t (SO) (50) 

will define SO galaxies in our model. Condensed objects, either having undergone few 
collisions, or which have collapsed too recently to form a disk, may be identified with 
damped Lya clouds, Lja forest or metal-line absorbers. 
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4.3. Conditions on the formation redshift of a galaxy 

From the modehng of sections §3.2 and 3.3 (equations (|30|) and (|33|)), we now rewrite 



conditions (^)-(|50|) as conditions on the formation redshift of galaxies. 



4.3.1. Field galaxies 
In the field limit, conditions (P5|)-(|5D|) become: 



1 + Zni<l + z{pi{z) (Spiral), (51) 

1 + Zni>l + ziii{z) (Elliptical) (52) 

and 

1 + 4,(z) <l + zni<l + zUz) (SO), (53) 

where zf^j and z^n are limiting redshifts obtained by setting A equal to A*^* and A^^', 

respectively, in equation (|30|). At 2; = (the present epoch), and provided Zni ^ 0, we find 
from (H) that 

l + ^4.(^ = 0)-(^)''' (54) 



and 



l + zUz = Q)^{^f\ (55) 



A. 

In general, z^^ and z^n are functions of the redshift z. 



4.3.2. Cluster galaxies 

Similar conditions on the redshift of formation can be obtained in the cluster limit from 
and (^)- (|50|) , provided that these conditions also hold in clusters. Those galaxies with 
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very large radii have stronger collisions and become ellipticals. Their redshift of formation, 
which at a given mass specifies their radius (pSj), is constrained by: 

Zni < zlii{M , p) (Cluster elliptical). (56) 

Spiral galaxies must be formed at redshift 

zlp^{M,p) < Zni (Cluster spiral) (57) 

and SO galaxies have a redshift of formation satisfying 

4;(M, p) < Zni < zlp^{M, p) (Cluster SO). (58) 

It is to be noticed that, contrary to the conditions in the field, cluster ellipticals form 
at a somewhat lower redshift than spirals. This is due to the competition between two 
effects: galaxies that formed earlier will have experienced more collisions by today, but 
recently formed galaxies will have larger radii and the efficiency of collisions will be higher. 
In clusters, the latter effect predominates as discussed in section §3.5: at a given mass, 
recently formed cluster galaxies have a higher collision factor than older ones. The limiting 
redshifts z^^j and z^n depend on the redshift z through (p3|). At the present time (z = 0), 
and provided ^ 0, we have from equation (^), with A'^'"''*'^'' equal to Aj^* and A^^^ 
respectively: 



Note that Af" and A^^* in (|5^) and (|60|) are not necessarily the same thresholds as in the 
field. 
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4.4. Summary 

We have presented a model for galaxy collisions. Using a phenomenological definition of 
galaxy morphological types, we have obtained, for each type in the limit z = 0, a condition 
on the redshift of formation of a galaxy: equations (|5lD -(^3D for field galaxies, (|56|)-(|58D for 



cluster galaxies. These conditions depend on the thresholds Afj^ and A^^*, which, as well as 
^^t, are free parameters of the model and scale all of the model predictions. 



5. Determination of 

The conditions derived above on formation redshifts are all independent of the initial 
spectrum of fluctuations, which need not be specified 0. Once an initial fluctuation spectrum 
is chosen, fluctuations at a given mass scale M are characterized by their effective height 
z/(M) in the matter density field linearly extrapolated until the present epoch: 

where Sc is the threshold on the linear density contrast required in structure formation 
theories, a the linear variance of mass fluctuations on scale M - equation and Zni the 
redshift at which fluctuations on scale M virialize. The mass function (0) can then be 
expressed in terms of the height u, given by (0), as: 

r/(M) = _ /2 A) cilna(M) r _ 1 . 



There is a consistency relation in order for the cooling constraint (^) that prevails 



at M ^ to provide the correct galaxy luminosity function ^{L). The local luminosity 



^We have introduced the CDM spectrum in §3.1 for the sole purpose of studying the 
variation of the mean galaxy mass M with redshift. 
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function of the Stromlo-APM redshift survey (Loveday et al. 1992) is well fitted by a 
Schechter function: 

$(L) = $r£)-exp(-L/L.)^, (63) 

where a ~ 0.97 and ~ 1.75x 10^^ Mpc^'^ (for h = 0.5). Assuming no luminosity evolution, 
we require that ( |63D and ( [6^ ) match on the scale M = M^, that is $(L)(iL|^, = r]{M)dM\^,. 
We get: 



p. 



z/* exp( 



(64) 



2 ' exp(l) V 2 po In (t(M) 
where z/^, = Sc{l + z*)/o"(M,,) and we have assumed a constant mass-to-light ratio. Solving 
equation (0), we find z/* ~ 2.8. Using CDM initial conditions (|TBp with a bias parameter 
6 = 1.67 ((78 ~ 0.6) required to reproduce the abundance of the observed X-ray clusters at 
present in an Qq = 1 universe (White, Efstathiou & Frenk 1993), we have a{M^) ^ 3.3, 
from which we infer 

(1 + 2*) ^5.5. (65) 

Once we know z*, the numerical values of the parameters of our model can be fixed. If we 
assume spherical symmetry, the radius of an = 10^^ Mq galaxy formed at redshift 
is K = (3M,/47rpol78(l + z^fy/^ ^ 50 kpc. Its velocity dispersion is ^ 0.5GM^/R^, 
implying K ~ 200 km/s. We can then compute the values of A* and Ac (equations (0) 
and (0)) provided we fix the quantities Vq and epo/M^,. The mean number density of 
galaxies today, epo/M^, is inferred from the local luminosity function of the Stromlo-APM 
survey, given by Loveday et al. (1992): epo/M^ ^ $* ^ 1.75 x 10"^ Mpc"^. Taking 
fo = 200 km/s, we find A* = 1.75 x 10~^ and A^ = 3.5 x 10^^. The numerical values of the 
parameters of the model are summarized in table 3. The only remaining quantities to be 
determined in order to fully normalize our model are A^^* and A^". 



- 25 - 



6. Determination of A^^* and A^" 

The quantities A^^* and A^" are fixed by requiring that the model produces the 
observed fractions of morphological types in the field today. These fractions can be 
evaluated in the following manner. Our phenomenological definition of morphological types 
implies definite conditions on the redshift of formation of galaxies. For example, field 
ellipticals are required to form at a redshift Zni > z^i (equation (|5^) ) in order to have 



experienced substantial energy exchange during collisions. The number density of objects 
of mass M which were already non-linear, with a mass between M/ A and M, at a redshift 
z' such that Zni > z' > z is given (Lacey & Cole 1993) by: 

n(M z z') - dlna~\M)6,{l + z) , 1/ ^.(1 + 



where 



and 



^ = ^-(^^ - ^) (67) 



2 

erfc(x) = —= / exp{—u^)du. (68) 



The erfc term represents the probability that an object present at z' still exists at redshift 
z. The present day fraction of ellipticals is thus given for z' = z^n by: 



/jjy™; r]{M, 0, ziii)dM 
i-,(M,0,OMA# (FWd elliptic* today). (69) 



The denominator of equation ([69| ) is the number density of all condensed objects with mass 
between Minf and Mgup, while the numerator counts only those objects that were formed 
before or at z*^" and that survived until the present epoch. Similarly, the fractions of spiral 
and SO galaxies are given by 



/it; [v {M,0,0)-v{M,0,z{ ^,)]dM 

jl^::;viM,o,o)dM 



^^vi = mIJ, ^ ^^1/^' (Field spirals today), (70) 
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and 

JM'-^[r]{M, 0, z{j - ri{M, 0, zfj]dM 

^ = -^^w ^p^^ ^"^^ (Field SOs today). (71) 

/£7^(M,0,0)dM 

Equations (|69|) -(|7lD depend on A^^* and A^" through the redshift cuts zf^j and zf^^. The 
observed fractions of morphological populations in the field today are ~ 65% spirals, ~ 10% 
ellipticals, and ~ 25 % SOs (Dressier 1980; Postman & Geller 1984). For the purpose of 
evaluating (|69D-([7TD, we use A = 2 in equation ( |67D and Mj„/ = 10^° Mq. Mgup is given by 
the cooling constraint (ESf). The CDM spectrum ([T6|) is used. We find in order to produce 



the correct abundance of morphological types that A^^* 0.003 and A^^' ^ 0.01. The 
collision factor for field ellipticals is thus typically larger by a factor at least 3 than the one 
for ellipticals. 



7. Redshift of formation of the morphological types in the field 

Equations (|5TD-(|53D give the typical redshift of formation of galaxies of different 
morphological types. Using the values derived previously for the thresholds on the collision 
factor, we find that typical field ellipticals form at Zni ^ 2.5 while typical field spirals form 
at somewhat smaller redshifts Zni ^ 1.5. SOs form, in our scenario, at intermediate redshifts. 
Figure 2 illustrates this prediction. The number density of the three galaxy types in the 
field today [z = 0) are plotted as a function of the redshift of formation of the galaxies. 
The results have been obtained from equations (|69l)-([n]) where the division by the total 
integrated galaxy number density has been omitted. The redshift cuts (pT|)-(pB|) have been 
used. The plot illustrates the main features of our phenomenological model: today's field 
ellipticals were formed at high redshift so as to experience efficient energy exchange through 
collisions. No galaxy forming at redshift lower than z ~ 2.5 will end up as an elliptical 
by today. Conversely, recently formed galaxies have suffered little tidal disturbance and 
constitute today's spirals. Finally, galaxies that were born at intermediate redshifts define 
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the present day SOs. 



8. Biasing between galaxy populations 

Each of the conditions (^)-(^) implies a precise value of the height z/(M) of 
fluctuations condensing out of the primordial density field into galaxies. For example, 
we infer from ( [5^ ) and ( [HT| ) that field ellipticals of mass M will condense out of linear 
fluctuations of height v such that 

where we define the threshold Vg^n for ellipticals. For field spirals of mass M, we have 

, , 1.68(1 + , , 

" ■= - MM) ■ 

where we define the threshold v^pi for spirals. For M^. galaxies, we have z/gj/ ~ 3 and z/^pj ~ 2. 

Equations ( [72|) and ( [73l) thus introduce a bias by effectively requiring that field ellipticals, 

spirals and SOs correspond to definite subsets of all mass condensations. By contrast with 

early theories of biased galaxy formation where v was free to be chosen, the formation 

threshold in our model results from the modeling of the physical processes involved during 

galaxy formation. 



9. The predicted morphology-density relation 

As early as 1958, Abell noted that elliptical galaxies are more frequently found in the 
cores of dense clusters (Abell 1958; Morgan 1961) while spiral galaxies predominate in the 
outer parts of clusters and make up almost 70% of field galaxies. Dressier (1980) first 
noticed the existence of a relation between galaxy morphological types and the (projected) 
local density in which galaxies are found. Postman & Geller (1984) and Giovanelli, 
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Haynes & Chincarini (1986) have extended Dressier 's work to assess the existence of a 
morphology- 3-D density relation in less dense environments such as for field galaxies and 
groups. Whether such a relation can be attributed to local properties of the environment, 
as claimed by Dressier, or can be related to initial conditions at birth (the 'nature versus 
nurture' controversy) has been an important issue assessed in more recent works (Whitmore 
& Gilmorc 1991; Whitmore, Gilmore & Jones 1993). Several effects have been proposed that 
could induce evolution of the morphology of galaxies. Spitzer & Baade (1951) pioneered 
such an approach by suggesting that spirals could evolve into SO galaxies by removal of their 
gas content. The dense cores of rich clusters provide an ideal environment: in these regions, 
ram-pressure stripping of the spiral interstellar medium by the intracluster gas is likely 
to occur (Gunn & Gott 1972). Direct mergers of disks (Mamon 1992) or tidal collisions 
(Spitzer & Baade 1951) are other likely effects. High resolution simulations suggest that 
compact sub-L^, SOs may form by the cumulative effect of tidal interactions on induced 
star formation and mass loss, so-called galaxy harassment (Moore et al. 1996). Competing 
theories have been developed by the advocates of the 'nature' hypothesis (Evrard et al. 
1990). These various approaches suffer however from a number of problems: while the 
Evrard et al. (1990) model fails to reproduce the observed morphology-density relation 
for SO and spiral galaxies, the result of removing the gas content of spiral galaxies by 
stripping would lead to SO-like galaxies, but not ellipticals, as stripping should not affect 
the stellar content of the stripped galaxy. Besides, Burstein (1979) observed that SOs 
seem to have thicker disks than spirals, which is difficult to reconcile with the hypothesis 
that spiral galaxies evolved into SO by gas stripping. The present model hes somewhat in 
between the 'nature' and 'nurture' hypothesis. It is of interest to test its predictions for the 
morphology-density relation. We undertake this task in this section. 

In order to evaluate the local morphology- density relation in our model, we need a 
continuous relation between the collision factor A and the density of the environment at 
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the present epoch. Scahngs of A with density have been obtained in the two hmiting cases 
of the field (equation (0), §3.2) and chister (equation (^7]), §3.4) environments. However, 
galaxies which are in a cluster now have not necessarily always been in the same cluster 
environment. In an = 1 CDM universe, clusters formed recently so that it seems fair 
to assume that energy exchange among colliding galaxies in a cluster supplement the early 
exchanges calculated for galaxies prior to cluster formation. In that case, the collision rate 
A for a given galaxy can be taken as 

A = A-^*''''^ + ^A'''"'**'''' (74) 

where A-^*^'*^ is given by equation (^) and A'^'"'**'^^ by equation (|33| ) using the scaling with 
the 3-D galaxy density n from equation (^Tj). In ([7^), we have introduced the dimensionless 



parameter 7 to take into account the fact that the thresholds A^^* and A^^', that define the 
various morphological types through the conditions (^8D -(pOD, may be different for cluster 
and field galaxies. According to ([7^), A-^**^'*^ dominates in the field, while the dominant 
contribution to the collision factor is A^'"**^** in high density environments. This defines 
a density n at which A*^'"***^^ starts to dominate over A-^**^'"^. Observationally (Postman 
& Geller 1984), the morphological populations are rather constant up to n ~ 1 Mpc~^ 
(/;,=0.5). Accordingly, we fix the value of 7 so that A^'"**^*" starts to dominate over A-^**^'*^ at 
n 1 Mpc~'^. Numerically, using ( |30D and (0), we find 7 ~ 10^'^. We have considered the 
case for an isothermal /5-model for the cluster gas density with (3 = 0.6 (Jones & Forman 
1984). The cluster velocity dispersion has been derived directly from the mass density 
profile using the virial theorem. Then, for n -C 1 Mpc~^, the proportions of the various 
morphological types are the same as in the field, while for n ^ 1 Mpc~^, they are derived 
essentially from A'^'"***^''. 

We can now evaluate the present day morphology-density relation, using equations 
(|69|)- ([7T|) where the redshift cuts z^u and Zgpi are now derived from ([7^) , with A = Afj^ and 
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A^^* respectively. Again, the CDM spectrum (|T6|) is used and is normalized to the present 
abundance of X-ray clusters (erg ~ 0.6). We use the mass function (|66D-(|67D with A = 2. 
Figure 3 shows the morphology- density relation obtained from our modeling (solid lines) 
compared to the observed relation (histograms; Postman & Geller 1984) at the present 
epoch. A typical l-cx error bar is shown for reference in the top panel (from Postman 
& Geller 1984). The model reproduces qualitatively the main features observed in the 
evolution with density of the proportions of morphological types: the spiral population 
decreases quite strongly in higher density environments, while this decrease is compensated 
by a corresponding increase of the elliptical and SO populations. The agreement is actually 
rather good quantitatively also, at least for low and intermediate densities. At high 
densities (n ~ 1000 Mpc~^), the model overestimates slightly the fraction of spirals, while 
the predicted fraction of ellipticals is somewhat lower than observed. 

Reproducing the observed morphology- density relation has required the adjustment of 
4 free parameters of our model: the two thresholds A^^* and A^"; the parameter 7 in (|74D; 
and the parameter (3 of the isothermal /3-model in ( ^7|) . The values of these parameters and 
the method used to determine them is summarized in table 4. The predictions of the model 
are quite sensitive to (3 and 7. For example, if we use the assumption that mass traces light 
in clusters, i.e. p ~ n, the model yields a decrease (an increase) of the spiral (elliptical) 
population at high density sharper than the ones observed. If we vary the value of 7, the 
density at which A'^'"'**'^^ starts to dominate over l^f^'^^'^ is different from the observed value 
n \ Mpc~^. It is however remarkable that the model simultaneously reproduces the 
morphology-density for the three Hubble types considered. Once the model reproduces the 
observed morphology- density relation today, quantitative predictions can then be made for 
the evolution with redshift of the various fractions of morphological types, in the field as 
well as in clusters. 
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10. Evolution of galaxy populations with redshift 
10.1. Redshift cuts 

10.1.1. Field Galaxies 



Conditions (^)-(50) on the collision factor A define the galaxy morphological types 
in our model. These conditions have been shown to be equivalent to conditions on the 
redshifts of formation of galaxies. The redshifts cuts on the redshift of formation, Zeii and 
Zspi, that delineate spirals from SO and from ellipticals are function of the redshift, since 
the collision factor itself depends on the redshift. Our model has been normalized so as to 
reproduce the observed fractions of the different galaxy populations today. It is a natural 
development of the model to predict the evolution of these populations with redshift. 

Equating A in equation (^) to the thresholds Af^ and A^^*, whose values have been 
determined previously, defines the two redshift cuts Zeu{z) and Zspi{z) at any redshift z. 
On figure 4, Zeii and Zspi are plotted (thin solid lines) as a function of redshift, for field 
galaxies. The field case is obtained by inserting the average galaxy density of the universe 
into equation ( |7^ ) so that A^'"***^^ is always small compared to A'^*^''^, whatever the redshift. 
The functions z^u^z) and Zspi{z) define three regions in the {z,Zni) space corresponding to 
specific ranges of values of A, as indicated on the figure. At a given redshift z, galaxies 
that went nonlinear at Zni > Zeu{z) are identified as ellipticals, those with Zni < Zspi{z) 
are spirals, and those whose redshift of formation lies in between the two redshift cuts are 
SOs. The thick solid arrow shows how a galaxy that formed at a given redshift Zni (here 
Zni = 3.1) would be identified at different redshifts: down to z ~ 2.7, this galaxy is identified 
as a spiral, as long as energy exchanges induced by tidal collisions are not high enough to 
meet our morphology criterion for SOs. From z ^ 2.7 to 2; ^ 1.8, the galaxy is seen as an 
SO, before becoming an elliptical when the cumulative effect of collisions since formation is 
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important enough to have induced substantial energy exchange. In our scenario, galaxies 
thus evolve along the Hubble sequence as time goes on. 

10.1.2. Cluster galaxies 

For cluster galaxies, the redshift cuts Zeii{z) and Zspi{z) are given by solving equation 
([7^) for A = Afil and A^^* respectively, this time with the average galaxy density of 
the universe replaced by the typical local galaxy density n of the cluster one wants to 
model. Figure 5 shows Zeu{z) for n = 500 Mpc~^ compared to the field case. For a given 
n, there is a redshift at which A-^*^''^ starts to dominate over 

^cluster _ At high redshift, 
/S^f^eid fiominates, as expected since clusters formed quite recently in a hierarchical universe. 
At low redshift, A^'"''*^'' dominates and Zeii{z) is different from the field values. In figure 
5, A'^'"**'^'' dominates over A-^*^''^ at z < 1.5. By construction, this redshift reflects the 
formation epoch of the cluster and is higher for higher values of n, since the cluster density 
reflects the density of the universe at the epoch of its formation. 

Another interesting feature in figure 5 is that the redshift cut Zeii is lower in a high 
density environment than in the field, implying that cluster ellipticals form at a lower 
redshift than in the field, as tidal collisions in clusters are more efficient than in the field at 
converting disks into ellipticals. 

10.2. Redshift evolution 

At any redshift, the knowledge of the cuts Zeu{z) and Zspi{z) allows one to count the 
number density of the different galaxy populations. As usual, we evaluate the fractions 
of the various populations from equations (|5PD-([7TD with ri{M,0, z^u) replaced, this time, 
by ri{M, z, Zeii), etc. All the results presented below assume the CDM spectrum ( |T6| ) 
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normalized to erg = 0.6. 

10.2.1. Definition of irregular galaxies 

So far, we have considered only three morphological types. However, observations 
at high redshift show the existence of a large population of so-called irregular galaxies 
(Brinchmann et al. 1997). These galaxies are morphologically perturbed, gas- rich, star 
forming galaxies with blue colors. There is a simple way to define 'irregular' galaxies in 
our model. Among all galaxies having experienced substantial energy exchanges during 
tidal collisions (essentially those that are ellipticals today), we can differentiate between 
those galaxies that have relaxed by the time characterized by redshift z and those that 
have not. The motivation for such a distinction is that unrelaxed galaxies may exhibit the 
disturbed morphology of observed irregulars, while relaxed galaxies are expected to show a 
smoother stellar distribution, as in ellipticals. Moreover, tidal collisions are likely to induce 
star formation in the interacting galaxies, and we can anticipate that the 'irregulars' of our 
model should be star-forming galaxies and have blue colors, in qualitative agreement with 
observations. However, part of the galactic gas is likely to be removed during collisions, 
and it is not clear whether 'irregulars' in our model will have retained enough gas to form 
stars. Clearly, modeling of the effect of coUisions on star formation and the physics of the 
gas is required. Such models are beyond the scope of the present paper and we shall restrict 
ourselves to the present definition. 

Our definition of 'irregular' galaxies thus amounts to introduction of a new redshift cut, 
Zirr- It is fixed by requiring that the elapsed time between the epoch of the last collision 
and the epoch characterized by z equals a typical time-scale for relaxation. Those galaxies 
that were formed at redshift higher than Zirr will have relaxed by redshift z and will be 
identified as ellipticals. In the opposite case, these galaxies will be 'irregulars'. In practice. 
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we assume that relaxation occurs on a time-scale treiax typically of the order of a dynamical 
time-scale tdyn, that is 

treiax ^ tdyn OC (X {1 + Znl)~^^'^ (75) 

where p is the density of the universe at the time of formation Zni- The epoch of the 
last collision is more problematic to evaluate. Actually, in our picture, they may not be 
any epoch of a last collision as collisions continuously occur over the lifetime of a galaxy. 
However, the bulk of collisions must occur at high redshift when the density of the universe 
is high. The redshift z^u (or a fraction of it) may be taken as the 'redshift of last collision', 
though more minor collisions may occur at lower redshift. There is some freedom in the 
choice of this epoch so that the definition of irregulars in our scenario is somewhat loose. 



10.2.2. Evolution of field galaxies 

Figure 6 shows the evolution of the population content for field galaxies. The main 
observed trends in the evolution of the galactic content of the universe are reproduced by 
our model: the fraction of field spirals decreases with increasing redshift as collisions in 
the past were more efficient at disrupting disk-forming galaxies. They make up 65% of the 
galaxy population today and only about 15% at redshift 3. The fraction of ellipticals rises 
smoothly from about 10% today to about 20% at redshift ~ 0.5, then starts decreasing 
while the fraction of 'irregulars' takes over rapidly and starts to dominate at redshift 2; ~ 1. 
We note that the substantial rise of the 'irregular' fraction at redshift z ~ 0.5 is, at least 
qualitatively, consistent with a recent analysis of the peculiar/irregular populations of the 
CFRS and LDSS redshift surveys (Brinchmann et al. 1997). 
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10.2.3. Evolution of cluster galaxies 

Figures 7 and 8 show the evolution of the galaxy population with increasing density of 
the environment. We take n = 10 Mpc~'^ in figure 7 and n = 500 Mpc~^ in figure 8. At 
high redshift, where A-^*'^^'^ dominates over A'^'"**^'" in equation (|7^ , the fractions of galaxy 
types are essentially the same as in the field. Going towards lower redshifts, the fraction of 
spiral galaxies increases less rapidly than in the field, peaks around z ~ 0.3 (figure 7) and 
steadily drops from 55% to 50% at the present time. This decrease is compensated by a 
corresponding increase of the SO and elliptical fractions, while the irregular population is 
quite unaffected. Thus, the cluster core is depleted from its spiral population as a direct 
consequence of the high density in the core, while it is populated by SO and ellipticals. In 
figure 8, similar behavior is noticed: the predicted fractions start to depart from the field 
fractions at redshift z ~ 1.5, as expected from inspection of figure 5, with the peak in the 
spiral fraction occurring at a somewhat larger redshift than in figure 7. The drop of the 
fraction of spirals is more pronounced and goes from 45% at redshift ~ 0.6 — 0.7 to 25% at 
the present time. It is temptative to associate the decrease of the fraction of spirals from 
intermediate redshift to the present with the well known Butcher-Oemler effect (Butcher & 
Oemler 1978; Dressier et al. 1994). Indeed, tidal collisions occurring shortly after cluster 
formation might induce star formation in spirals, that then would appear similar to the 
blue Butcher-Oemler galaxies, while exhaustion and/or stripping of the gas content during 
subsequent collisions would transform spirals into SOs and ellipticals. Clearly, inclusion of 
collision induced star formation and gas physics in our model is needed for any detailed 
investigation of its predictions with respect to the Butcher-Oemler effect. 



-se- 
ll. Conclusion 

We have argued that galaxy coUisions play a fundamental role in galaxy formation. We 
have built a model in which galaxy morphological types originate as a result of gravitational 
interactions with surrounding galaxies. Simple rules for energy exchange during collisions 
have been proposed that allow us to discriminate between different Hubble types: efficient 
collisions result in the disruption of disks and substantial star formation, leading to the 
formation of ellipticals; few or inefficient collisions in the past allowed a large gas reservoir 
to survive and form disks. Quantitative analysis of energy exchanges in the field and cluster 
environments have been presented. These analyses are based on the simulation results of 
Aguilar & White (1985). 

Assuming that galaxy formation proceeds in a = 1 Cold Dark Matter universe, our 
model both reproduces a number of observations and makes various predictions, among 
which are the redshifts of formation of the different Hubble types in the field. When 
the model is normalized to the present day abundance of X-ray clusters, we find that 
the amount of energy exchange needed to produce ellipticals implies that they formed by 
2; ^ 2.5 while spirals formed at z ^ 1.5. 

The model also provides a natural basis for biasing between field spirals and eUipticals 
without requiring an ad-hoc identification of morphological types to peaks of different height 
in the initial density field. Field elliptical galaxies are found to be more biased with respect 
to mass than spiral galaxies by a factor ~ 1.5. They preferentially form out of ~ 3a peaks, 
whereas spiral galaxies condense out of ~ 2a peaks. 

Our formalism allows us to study galactic evolution in clusters. With the same collision 
rules as in the field, the model satisfactorily reproduces the morphology- density relation 
which spans a range of densities from 10~^ to 10'' galaxies per Mpc^. 
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Finally, the predictive power of the model is exploited to predict the evolution of 
galaxy populations with redshift. The predicted trends are in good qualitative agreement 
with observations, both in the field and in clusters. Our modeling of collisions in clusters 
naturally gives rise to a Butcher-Oemler-like effect, observed in the predicted spiral 
population. However, inclusion in our model of gas physics and collision induced star 
formation is required for any detailed quantitative comparison of the predictions of our 
model to observations. Such a comparison will be undertaken in a forthcoming paper. 

We thank G. Mamon, A. Blanchard, P. Ferreira, B. Moore, S. Zepf and R. Wyse 
for fruitful discussions or comments regarding the results presented in this paper. 
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been supported in part by a grant from the NSF. 



-38- 

REFERENCES 

Abell, G. O. 1958, Astr. A&AS, 50, 241 

Aguilar, L. A., & White, S. D. M. 1985, ApJ, 295, 374 

Bahcall, N. A. 1979, ApJ, 232, 689 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 
Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 

Bingelli, B. 1989, in Large scale structure and motions in the universe, ed. M. Mezetti 
(Dordrecht), 47 

Blumenthal, G. R., Faber, S. M., Primack, J. R., & Rees, M. J. 1984, Nature, 311, 517 

Brinchmann, J., et al. 1997, preprint 

Burstein, D. 1979, ApJ, 234, 829 

Butcher, H., & Oemler, A., JR. 1978, ApJ, 226, 559 

Cavahere, A., & Fusco-Femiano, R. 1976, A&A, 49, 137 

Clements, D. L., Sutherland, W. J., McMahon, R. G., & Saunders, W. 1996, MNRAS, 279, 
477 

Cole, S., Aragon-Salamanca, A., Frenk, C. S., Navarro, J. F., & Zepf, S. E. 1994, MNRAS, 
271, 781 

Dekel, A., & Silk, J. 1986, ApJ, 303, 39 
Doroshkevich, A. G. 1970, Astrofisika, 6, 581 
Dressier, A. 1980, ApJ, 236, 351 



Dressier, A., Oemler, A., JR., Butcher, H. R., & Gunn, J. E. 1994, ApJ, 430, 107 
Evrard, A. E. 1989, ApJ, 341,26 

Evrard, A. E., Silk, J., & Szalay, A. S. 1990, ApJ, 365, 13 

Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 

Giovanelli, R., Haynes, M. R, & Chincarini, G. L. 1986, ApJ, 300, 77 

Gunn, J. E., & Gott, J. R. 1972, ApJ, 176, 1 

Jones, C., & Forman, W. 1984, ApJ, 276, 38 

Lacey, C. G., & Cole, S. 1993, MNRAS, 262, 627 

Lacey, C. G., & Fall, S. M. 1983, MNRAS, 204, 791 

Leech, K. J., Rowan-Robinson, M., Lawrence, A., & Hughes, J. D. 1994, MNRAS, 267, 675 
Lemaitre, G. 1933, Ann. Soc. Sci. Bruxelles, A53, 51 

Loveday, J., Peterson, B. A., Efstathiou, G., & Maddox, S. J. 1992, ApJ, 390, L338 
Mamon, G., A. 1992, ApJ, 401, L3 

Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A., JR. 1996, Nature, 379, 613 
Morgan, W. W. 1961, Proc. Nat. Acad. Sci, 47, 905 
Narayan, R., & White, S. D. M. 1988, MNRAS, 231,97 

Peebles, P. J. E. 1980, The large-scale structure of the universe (Princeton Univ. Press) 
Postman, M., & Geller, M. J. 1984, ApJ, 281, 95 
Press, W. H., & Schechter, P. L. 1974, ApJ, 187, 425 



-40- 

Rees, M. J., & Ostriker, J. P. 1977, MNRAS, 179, 451 
Richstone, D. O. 1975, ApJ, 200,535 
Richstone, D. O. 1976, ApJ, 204,642 

Rocca-Volmerange, B., & Schaeffer, R. 1990, A&A, 233, 427 

Schaeffer, R., & Silk, J. 1985, ApJ, 292, 319 

Schweizer, F. 1996, AJ, 111, 109 

Schweizer, F., & Seitzer, R 1992, AJ, 104, 1039 

Sellwood, J. A., & Carlberg, R. G. 1984, ApJ, 282, 61 

Silk, J. 1977, ApJ, 211, 638 

Spitzer, L. 1958, ApJ, 127, 17 

Spitzer, L., & Baade, W. 1951, ApJ, 113, 413 

Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 

Toomre, A. 1974, in The formation and dynamics of galaxies, ed. D. Reidel (Dordrecht), 
347 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 
White, S. D. M. 1984, ApJ, 286, 34 

White, S. D. M., Efstathiou, G., & Frenk, C. S. 1993, MNRAS, 262, 1023 
Whitmore, B. C., & Gilmore, D. M. 1991, ApJ, 367, 64 
Whitmore, B. C., Gilmore, D. M., & Jones, C. 1993, ApJ, 407, 489 



-41 - 



Figures captions 

Fig. 1. — Average galaxy mass M (in Mq) as a function of redshift, with and without the 
coohng constraint (23). When coohng is included, M is relatively constant up to 2; ~ 4. 

Fig. 2. — Number density of present day field galaxy populations as a function of the redshift 
of formation of galaxies. This graph shows the contribution of the three morphological types 
defined in our model to the present population of galaxies. Galaxies formed at redshift 
z > z(ii become ellipticals by today. Galaxies formed at z^^^ < z < z^n are present day 
SOs, while galaxies forming at z < z{p^ arc today's spirals. The graph illustrates how 
morphological populations are defined and counted from a combined use of the mass function 
(66) and the predicted redshift cuts z^n and zf^j 

Fig. 3. — Predicted local morphology-density relation (sohd hues) compared to the data 
(histograms) of Postman & Geller (1984). A typical one standard deviation is shown on 
the top panel. The morphology-density relation is computed in our model by assuming that 
energy exchange in a cluster A'^'"**'^'' supplement the early exchanges calculated for galaxies 
prior to cluster formation A-^*^'*^. Scahng of A'^'"**^'' with the local galaxy density n is obtained 
by assuming that the cluster density distribution traces the gas density distribution. An 
isothermal /3-model is assumed for the gas distribution. 

Fig. 4. — Redshift cuts for the formation of Hubble types, z^^ and z{pj^ (thin solid lines), are 
plotted as a function of redshift in the field, z^n and ^f^j define three regions in the {zni, z) 
space, corresponding to three ranges of values of the coUision factor A: eUipticals lie in the 
upper region. Spirals in the lower region, while SOs lie in the intermediate region. The solid 
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thick arrow shows how the morphological type of a galaxy that was born at redshift 3.1 
evolves with redshift, under the effect of coUisions. 

Fig. 5. — The redshift cut z^i in the cluster environment (dashed hne) is compared to the 
same cut in the field (solid line). At high redshift, the collision factor (74) is given by the 
field limit A^^'^^'^, while at lower redshift A'^'"'**'^'' starts to dominate over A^^'^^'^ and significant 
discrepancy with the field appears. Ellipticals in a cluster environment are able to form at 
lower redshift than in the field due to the higher efficiency of collisions in clusters. 

Fig. 6. — Predicted fractions of galaxy populations in the field as a function of redshift. A 
new morphological type has been introduced and labeled as 'irregulars' (see text). 

Fig. 7. — Predicted fractions of galaxy populations for n = 10 Mpc~^ as a function of 
redshift. At high redshift, predictions are the same as in the field, whereas, at low redshift, 
significant differences appear due to the different physics of collisions in denser environments. 

Fig. 8. — Same as figure 7 for n = 500 Mpc~^. Increasing n accentuates the features of the 
evolution of the galaxy population at low redshift. 
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Table 1: Name and meaning of parameters used 



Name 


Meaning 


R 


Radius of the test galaxy 


M 


Mass of the test galaxy 


V 


Internal stellar velocity dispersion of the test galaxy 


V 


Relative velocity of encounters 


e 


Mass fraction of the universe in galaxies 


P 


Average density of the universe 


f 


Normalization of the collision efficiency of AG85 function Je 


A 


Rate of change of binding energy due to collisions 



Table 2: Rate of change of binding energy 



Environment 


Equat 


ion Used 


Underlying Assumption 


Cluster (v > V) 


(0) 


Impulse Approximation 
Straight line trajectories 


Field (v^V) ( 


11 


) Gravitational focusing 
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Table 3: Numerical values of the model parameters 



Parameter 


Numerical Value 


Units 




1 




Ho 


50 


km/s/Mpc 




R* 


50 


kpc 


K 


200 


km/s 




4.5 




epo/M^ 


1.75 X 10-3 


Mpc-3 




200 


km/s 


f 


30 




A, 


1.75 X 10-^ 




Ac 


3.5 X 10-^ 





Table 4: Tuned parameters for the morphology- density relation 



Parameter 


Numerical Value 


Fixed By 




0.6 


Isothermal /3-model 


7 


10-3 


Galaxy populations at n = 1 Mpc-^ 


Ar 


0.003 


Galaxy populations 




0.01 


in the field today 



T 1 1 1 i 1 1 □ 



10 



14 _ 



Average mass [Equation (12)] 



10 



13 



No Cooling 



1012 



10^1 



With Cooling 



10 



10 



J I I I L 







1 2 3 4 5 6 7 8 

Redshift z 




12 3 4 



Redshift of Formation 



0.8 
0.6 
0.4 
0.2 

o 

? 0.8 

CO 

c£ 0.6 
•2 0.4 

2 0-2 
ex 

o 

^ 

0.8 
0.6 




Ellipticals 



0.4 



0.2 




J I I I I I I I I 





0.01 0.1 



1 10 100 1000 

n [Mpc"3] 



